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Recent experiments with a rotating probe immersed in weakly fluidized granular materials show a 
complex behavior on a wide range of timescales. Viscous-like relaxation at high frequency is accom¬ 
panied by an almost harmonic dynamical trapping at intermediate times, with possibly anomalous 
long time behavior in the form of super-diffusion. Inspired by the Itinerant Oscillator model for dif¬ 
fusion in molecular liquids, and other models with coupled thermostats acting at different timescales, 
here we discuss a new model able to account for fast viscous relaxation, dynamical trapping and 
super-diffusion at long times. The main difference with respect to liquids, is a non-negligible cage 
inertia for the surrounding (granular) fluid, which allows it to sustain a slow but persistent mo¬ 
tion for long times. The computed velocity power density spectra and mean-squared displacement 
qualitatively reproduce the experimental findings. We also discuss the linear response to external 
perturbations and the tail of the distribution of persistency time, which is associated with superdif¬ 
fusion, and whose cut-off time is determined by cage inertia. 

PACS numbers: 45.70.-n,05.40.-a,47.57.Gc 


I. INTRODUCTION. 


Granular materials share analogies with condensed 
“molecular” matter, but often escape its well-established 
theoretical approaches 0 0- Equilibrium statistical 
physics may suggest only very approximate ideas about 
the qualitative behavior of granular media in strongly 
fluidized steady states and dramatically fails in the ex¬ 
treme case of static or quasi-static regimes. Continuum 
descriptions for dense flows lack first-principle constitu¬ 
tive relations [3, jJ, while more refined kinetic theories 
(e.g. mode-coupling) must be carefully adapted to take 
into account some fundamental peculiarities, such as the 
breaking of time-reversal invariance and the relevance of 
the inertia of the medium jHj. The building up of granu¬ 
lar hydrodynamics from “microscopic” models (ie. where 
all grains are described) is a promising program, but for 
the moment its success is limited to dilute regimes and 
rests upon the (uncertain) separation between fast and 
slow scales @-0- 

Our understanding of the liquid state of granular mat¬ 
ter, being in the middle between two opposite worlds (the 
very fast “granular gases” and the very slow “granular 
glasses”), is even more incomplete and may benefit from 
simplified effective models. An important insight is pro¬ 
vided by experiments, where the “liquid” state is realized 
through the application of some mild shaking leading to 
a slowly mixing flow with strong correlations and long 
but finite relaxation times M- When the longest re¬ 
laxation time overcomes the experimental times, one may 
say to have reached a transition point , entering into a sort 
of - empirically defined - solid or glassy state [12ijl4j . We 
do not intend to directly address such a transition: how¬ 
ever some of our results, in the following, concern also 
this delicate point. 

The present paper aims at discussing a simplified lin¬ 
ear model which is able to reproduce some noticeable 


phenomena observed in a granular liquid state 11]. In 
particular our ambition is to propose a minimal model 
which exhibits a transient cage effect and super-diffusion 
at later times. Caging is a common hallmark of diffusion 
in dense liquids [15( and it is usually found in granular 
systems at large packing fractions 0 ; 16|, |l3] ■ Superdif¬ 
fusion is much less common in liquids, it seems rather 
a peculiar effect of granular systems [ 13 , 0: however 
it is rarely seen and hardly explained: below the jam¬ 
ming transition it has been observed in fl4| . above such 
a transition it was seen in fl3| where it was imputed 
to “zero”-modes of the host fluid, or in |l8| where the 
mechanism of Taylor dispersion was involved, or in jl9j 
explained by a turbulence-like cascade effect. A universal 
scenario for anomalous diffusion is lacking [20|, but cer¬ 
tainly it is the signal of an enduring memory. A family 
of phenomenological models for anomalous diffusion in¬ 
cludes fractional Fokker-Planck equations 21], where an 
immediate physical interpretation is not always at hand. 
The observations in EH were better explained through 
a phenomenological continuous time random walk model 
for the velocity [20 ] . with a power-law-decaying distri¬ 
bution of persistency times which was confirmed by ex¬ 
perimental measurements. Such a model however could 
not explain the cage effect (which is a sub-diffusive be¬ 
havior at earlier times) and was, therefore, adopted to 
match only partially the experimental results, in partic¬ 
ular the slow time-scales. The model presented here, on 
the contrary, aims at offering a unifying picture for the 
two phenomena, and highlights the essential role of the 
“cage inertia”, which is the origin of long-time memory. 

The organisation of the paper is as follows. In Sec¬ 
tion II we summarize the results of a recent experiment, 
carried on by one of the authors, where time correla¬ 
tions and mean squared displacement of a probe were 
measured, interpreted under the light of a first simpli¬ 
fied model. In Section III we propose the new Itinerant 
Oscillator model with cage inertia, with a discussion of 
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its motivation. In Section IV we report the main for¬ 
mula for static quantities in the steady state. In Section 
V, we discuss two-time quantities in the steady (time- 
translantion invariant) state, including the velocity spec¬ 
trum, the mean squared displacement, linear response 
and, computed only numerically, the distribution of per¬ 
sistency times. Conclusions and perspectives are drawn 
in Section VI. 


II. A RECENT EXPERIMENT ON GRANULAR 
MESORHEOLOGY 

The granular liquid state is characterized by the emer¬ 
gence of many time-scales, associated with the complex 
and collective relaxation behavior of grains. A window 
into those time-scales may be open by studying the dy¬ 
namics of a diffusing impurity, both in experiments and 
in simulations [III [HI • A recent experimental study f 1 1 ] 
has offered a new picture in a wide range of time-scales, 
from 10 -3 s up to 10 3 s and more, revealing a non 
trivial scenario. In the experimental setup, sketched in 
Fig.[U\ the “impurity” was constituted by an immeresed 
blade (with momentum of inertia I) who could rotate 
around a fixed vertical axis under the kicks from the 
grain of a vibrofluidized granular medium. The dynam¬ 
ics of the angular velocity u(t) of the blade and its ab¬ 
solute angular position 6(t) = f 0 dscu(s), was studied in 
different regimes of density and intensity of vibration. 
In Fig. [TJ3, the velocity power density spectrum (vpds) 

S(f ) = 2ntloT l/o r ° T ^{t)e i{2nf)t dt\ 2 is presented and its 
salient features are highlighted in two opposite limits, 
which are the gas and the cold liquid. We remind that 
the vpds is the Fourier transform of the velocity autocor¬ 
relation function (vacf) and that its / —?> 0 + limit is the 
self-diffusion coefficient, i.e. D x = 7rlimy_ >0 + <?(/)• We 
also recall that relations exist, under certain approxima¬ 
tions, between the vpds and the intermediate scattering 
function which - in liquids - is typically accessed through 
neutron scattering experiments [241 ] . 

In the gas limit, when the packing fraction is low and 
the average energy per grain is high, the probe veloc¬ 
ity autocorrelation (vacf) is close to a simple exponen¬ 
tial decay ~ e~ t / Taas , ruled by a single relaxation time 
T gas : in this limit the vpds takes the form of a Lorentzian 

S(f) « ^/I 1 + ( 2 tt/// 7 ) 2 ] with 7 = I/r gas . 

In the cold liquid limit, when the packing fraction is 
high (larger than 30 — 35%) and the average energy per 
grain is small (but still in an “ergodic” phase), the ob¬ 
served vpds strongly deviates from the Lorentzian. Ig¬ 
noring a mechanical resonance due to the mounting plate 
at ~ 10 2 Hz, it displays four different regions: at high fre¬ 
quency (region IV) it decays with a negative power law 
equal or smaller than 2; in region III it shows a smooth 
parabolic maximum (centered near ~ 10 Hz, reminiscent 
of a harmonic confinement (“cage”); in region II it stabi¬ 
lizes on a short plateau, which suggests a loss of memory 
(as in the plateau of the Lorentzian which marks the on- 
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FIG. 1. A: sketch of the experiment reported in 0. B: 
experimental data of the vpds for the gas case and the “cold 
liquid” case, together with predictions (dashed lines) from the 
incomplete model, Eq. ©• C: experimental data of the msd 
for both cases, together with dashed lines useful as guides for 
the eye. 


set of normal diffusion); finally region I, perhaps the most 
surprising one, shows an increasing S(f) for / —> 0 + , sig¬ 
naling a problem with the finiteness of the self-diffusion 
coefficient D^. A few longer experiments (12 hours) were 
conducted, showing a slow crossover toward a new higher 
plateau at very low frequencies. The study of the mean 
squared displacement (msd), see Fig. [Tp confirmed that 
the four regions of the cold liquid case correspond, re¬ 
spectively, to short-time ballistic (free) motion (IV), dy¬ 
namical arrest due to caging (III), later relaxation of the 
cage (II) and “final” superdiffusive behavior (I). 

In [llj a first model was proposed to account for regions 
II-III-IV of the vpds in the cold liquid limit. The model, 
which was inspired by a model for diffusion in cold liquids 
called Itinerant Oscillator model [25140 ]. also recently 
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used to describe microrheology in living matter j28j . de¬ 
scribes the evolution in time of the angular velocity of 
the probe w(f) and its angular position 9(t) = f* dsoj(s), 
according to the following stochastic equations of motion: 


Iio{t) = -yw(t) - k[9(t) - 9 0 (t)] + (la) 

6(t) = f w(s) (lb) 

Jo 

e 0 (t) = 2D 0 f rio(s) (lc) 

Jo 

(Id) 


where rj(t) and rjo(t) are independent white normal Gaus¬ 
sian noises (unitary variance). The model represents the 
diffusion of a particle in a harmonic potential with “stiff¬ 
ness” k and unfixed minimum located at 9o(t), under the 
effect of a thermal bath at temperature T and relaxation 
time I/'y. The harmonic potential, representing the cage 
created by the confining effect of the dense granular host 
fluid, is not fixed but moves, as $o(t) behaves as Brown¬ 
ian motion with diffusivity Dq. Motivation for this model 
is twofold: 1) it reproduces the main features of the vpds, 
i.e. short time fast relaxation (region IV), an elastic res¬ 
onance at intermediate times (region III) and a plateau 
revealing loss of memory at larger times (region II); 2) 
in the dilute limit (when k —> 0) it can be rigorously 
derived [29| , 3) at intermediate densities a series of stud¬ 
ies showed that memory effects (coming from correlated 
collisions) are well described by a similar coupling with 
an additional degree of freedom characterized by slower 
relaxation time-scales [23| . The vpds of the above model 
can be calculated and reads 


S(f) = - 


D 0 k 2 + 1 T(2nf) 2 

212 


7r 7 2 (27 t/) 2 + [k — J(27t/) 2 ] 


( 2 ) 


Two limiting cases are recovered: when k = 0, the 
Ornstein-Uhlenbeck process is obtained, with S(f) tak¬ 
ing the Lorentzian form mentioned before. When k > 0 
and Dq = 0, one has the Klein-Kramers process in a fixed 
harmonic potential, and S(f) —» 0 for / — > 0, expressing 
the absence of diffusion at large times: the cage does not 
move and fully confines the particle. Formula ([2]) fairly 
fits the experimental spectra (see dashed lines in Fig.[lj3) 
in regions II-IV, with k/I ~ (27r • 10) 2 ~ 4 • 10 3 Hz 2 see 
El- Reasonably, the “cage stiffness” decreases at in¬ 
creasing shaking intensity. It also decreases as the density 
is reduced, and abruptly goes to zero at packing fractions 
of the order of ~ 15%. The “cage diffusivity” D 0 rapidly 
increases with increasing T and with decreasing packing 
fraction (or number of particles). 

The main problem of model m is that in region I it 
always predicts a diffusive behavior, no super-diffusion is 
allowed. In El superdiffusion was put in strict relation 
with a broad distribution of “inversion”, also called “per¬ 
sistency” , times measured with the following recipe. The 
frequencies of regions III and IV were filtered out, by con¬ 
sidering a smoothed angular velocity oo s (t) = i uj(s)ds 


(t > Is), which displayed much smaller fluctuations 
than those of w(t), but with positive correlation for very 
long times (larger than 10 s). This behavior, incompat¬ 
ible with model ([!]) (whose relaxation times are much 
smaller), is likely to be due to a quasi-steady motion of a 
large part of the granular medium surrounding and there¬ 
fore dragging the probe. The large inertia is rensponsible 
for the observed long relaxation times. At high frequen¬ 
cies (part of region II, region III and IV) the contribution 
of cu s (t) is negligible, explaining the good agreement with 
model © in that part of the spectrum. 

In the following Section we propose an extension of 
this model, taking account the inertia of the surrounding 
medium, in order to reproduce the superdiffusive behav¬ 
ior. 


III. THE ITINERANT OSCILLATOR MODEL 
WITH CAGE INERTIA 

Motivated by the experimental measurement of long 
relaxation times, we introduce a new model to replace 
(or extend) that of Eq. Jl]): 

Iu(t) = - 7 oj(t) - k[9(t ) - 0 o (t)] + ^/2jTr](t) (3a) 

/ 0 w 0 (t) = - 70^0 W + k[9(t) - 9 0 (t )] + \j2'yoTorjo (i) 

(3b) 

6(t) = ( Lo(t')dt' (3c) 

Jo 

e 0 (t) = [ u>o(t')dt'. (3d) 

Jo 

In the above equations r](t) and rjo (t) are independent 
white gaussian noises with zero average and unitary vari¬ 
ance, namely ( rj(t)rj(t')) = S(t — t') and (: r] 0 (t)r](t')) = 
5(t - t'). 

In equations Q the angular velocity of the probe 
feels two different forces: of course they are both re¬ 
lated to collisions, but one part is without memory 
and is described by the Ornstein-Uhlenbeck contribu¬ 
tion — 7 ui + T](t), while the second part takes the form 
—k[9(t) — 9o(t)] and therefore depends upon the past his¬ 
tory of u(t) and u>o(t). The choice of a harmonic inter¬ 
action aims at simplifying the computations and can be 
justified by the small velocities of the blade with respect 
to that of the surrounding particles. The delaying force is 
modelled as a drag toward a reference point which slowly 
evolves in time, according to Eqs. m-m- The variable 
9o(t) should be viewed as a collective degree of freedom 
representing the preferential point of the blade with re¬ 
spect to some granular cage. The cage slowly changes its 
configuration and favours the blade’s drift at later times. 

In the previous, simpler, model, Eq. Q, the dynamics 
of 6o(t) was overdamped , as in a motion without inertia. 
The new model takes into account the crucial effect of 
cage inertia, through the introduction of a fourth degree 
of freedom ojq (t) = 9 q which evolves with Eq. © , as 
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well as the effect of the blade upon the granular material 
through the reciprocal elastic drag. This last ingredient 
is likely to be negligible, in view of the large value of the 
inertia Jo, but its inclusion is convenient for symmetry 
and physical consistency. We must emphasize that this 
model can be understood in the context of other mod¬ 
els, often referred to as “two temperature models” (see 
|30l (3l| and references therein), where a particle moves 
under the influence of different thermostats acting at dif¬ 
ferent timescales. Even if it is quite natural, in dissipative 
systems such as a granular fluid, to associate different 
temperatures to different timescales, in the following we 
show that the main subject of our study, which is super¬ 
diffusion, can be obtained even for equal temperatures. 



1 

T 

7 

k 

Io 

To 

70 

case A 

1 

1 

too 

0 

/ 

/ 

/ 

case B 

1 

1 

too 

5000 

0 

0.1 

1000 

case C 

1 

1 

too 

5000 

10 

0.1 

1000 

case D 

1 

1 

too 

5000 

10 3 

100 

10 

case E 

1 

1 

too 

5000 

10 4 

100 

10 

case F 

1 

1 

too 

5000 

10 4 

1 

10 


TABLE I. Values of the parameters for the cases illustrated 
in the paper. A: The Ornstein-Ulenbeck process. B: A case 
without cage inertia (overdamped cage). C, D, E: cases with 
small, medium or large cage inertia. F: equal to case E, but 
at thermodynamic equilibrium To = T. 


Introduction of cage inertia is the main novelty with 
respect to the original Itinerant Oscillator model, includ¬ 
ing the version discussed in [ill ], and certainly deserves 
some motivation. At low temperatures, activated pro¬ 
cesses (that is the possibility to “jump out” of the cage 
thanks to some thermal fluctuation) are negligible and 
the probe never really escapes from a cage: on the con¬ 
trary, it is the cage that slowly evolves and dictates the 
motion of the probe at large times. A displacement of 
the probe A9(t) between two instants separated by a 
large time t , therefore, is closely related to a displace¬ 
ment A9 0 (t) of the cage itself, i.e. of a large part of 
the surrounding granular medium. Then it is reasonable, 
when looking for an ingredient reproducing almost bal¬ 
listic superdiffusion, to inunagine that the cage is doing 
long ballistic drifts (at very small velocity), sustained 
by its large inertia. Discreteness and finiteness of the 
granular material, which in the experiments is made of a 
few thousands grains, makes random but persistent (also 
called “secular” [Hj]) drifts possible. For those reasons it 
seems reasonable to us the introduction of a large “cage 
inertia” Iq. The event of a large portion of solvent to 
drift in a particular direction for times larger than the 
cage relaxation is highly unlikely in a dense molecular 
liquid: this explains why one usually do not observe su¬ 
perdiffusion in ordinary liquids. It is out of the scope 
of the present paper to discuss the exact mechanism of 
formation of the large cage inertia Iq in a vibrated gran¬ 
ular medium, as well as the reason why it increases when 
the vibration amplitude is reduced [33J. The answers 
to such questions are postponed to future investigations: 
our main aim, here, is to convince the reader that the 
concept of “cage inertia” is useful for the description of 
low temperature granular liquids. 

The model is linear and it is useful to recast it into 
a more compact form, by defining the vector w(t ) = 
{z(t), w(i), wo(t)}, with z(t) = 9(t) — 9 0 (t). The model 
then takes the form 


( 4 ) 


with 


A = 


B = 


V = 





( 5 ) 

( 6 ) 

( 7 ) 


The analysis of the model will be discussed focusing 
on a few particular cases, inspired by the experimental 
results of [ll| , which are listed in Table Q] By fixing 7=1 
and T = 1, the units of moment of inertia and time are 
fixed. A comparison with experimental observations sug¬ 
gests that our arbitrary time unit is close to 1 real sec¬ 
ond. In Table [H case A is similar to a dilute experiment, 
while cases D or E are similar to a dense and cold one. 
Case B is an example of cage without inertia (Iq = 0), 
which cannot display superdiffusion. Case C is similar 
to that, as Iq is small. Finally Case F has still a large 
cage inertia, but has the peculiarity to be at thermody¬ 
namic equilibrium, i.e. Tq = T . Nevertheless, we will 
see that this ingredient is relevant only in the study of 
Fluctuation-Dissipation relation and has no consequences 
for the presence of cages or superdiffusion, as it does not 
crucially affect the timescales of relaxation. 


IV. STATICS: ONE TIME QUANTITIES IN THE 
STATIONARY STATE 


Eqs. 0 constitute a linear system, which is solved by 
Gaussian multivariate distributions. We assume that the 
values of the parameter of the model are such that only 
eigenvalues with positive real part appear and a steady 
state can be reached. In such a steady state, the one time 
distribution function takes the form 


P(z,uj,ujo ) 


1 

yj (2n) 3 det(a) 


Exp 



w = —Aw + Bfj(t), 


( 8 ) 
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where (3 = a 1 and a is the covariance matrix, which we 
write in the form 


&ZZ &ZUJ &ZUJQ 
® ® ZUJ &UJUJ &OJOJO 

y &zoj a &OJOJQ ^ujqujq 

The covariance matrix satisfies the equation [341 
Ao + crA T = BB, 
which is rewritten as 



& ZUJ — O ZUJ 0 

I {^UJUJq G UJUJ ) + k&ZZ H"~ 70 ZUJ — 0 
-^0 ^UJqUJq &UJUJQ ) k&ZZ H - ^O&ZUJo — 0 

I{ka zu] + ^(Ju>u) = T'y 

Io(.10&UoUJO o) ^070 

k@Z0J0 4 " I^UJUJQ _j_ 'YO&UJUJO fo&ZLJ _ „ 

7 To ~ • 


The solution reads 


&UJUJ 




Ozz 


&UJUJQ 
& ZUJ 


^+7 (T-T,)A 

Io d 

^ N 4l3 


( 9 ) 


( 10 ) 


(lla) 

(lib) 

(11c) 

(lid) 

(lie) 

(Ilf) 


(12a) 

(12b) 

(12c) 

(12d) 

(12e) 


With Ax = k( lo I + 7 / 0 ), A 2 = 7o(/- fc (7 + 7c)+7So/ + 7 3 /o) i 
A3 = 770 k(I - lo), A 4 = 770 ( 70 1 + 7 / 0 ) and d = (7 + 
7o) ( 70 1 2 k + 77 qI + 7^0 + 7 2 7o Io) ■ 

We can recover two well known physical limits. In the 
decoupling limit k —> 0, one has 


_ T 
&ujuj — ~Y 

- 

°7’0“0 — T 
JO 

CT zz = OO 

^UJUJq 0 

& ZUJ & ZUJ Q 


T-Tp 
7 + 70 ’ 


(13a) 

(13b) 

(13c) 

(13d) 

(13e) 


It is important to note that this limit is singular (since the 
case k = 0, corresponding to two independents Ornstein- 
Uhlenbeck processes, is non-stationary for z): this ex¬ 
plains why a ZUJ = a ZUJo ± 0. 

For T = Tq, that is at thermodynamic equilibrium, one 



FIG. 2. We show how an (for some i, i.e. uj, ujq and z, 
properly normalized) depend on k or Jo, while all the other 
parameters are fixed according to case C of Table Q] On the 
top panel, we see the important effect of the coupling quan¬ 
tity k on both self-correlations and cross-correlations. On the 
bottom one, we notice that the effect of cage inertia Io on the 
coupling (e.g. cr„„ 0 , green curve), is larger for intermediate 
values: at smaller values cr„ 0 „ 0 is close to its “uncoupled” 
value, while it deviates from it at large values; the opposite 
happens to a wu . 


has 


T 

(14a) 

&ujuj — y 

T 

T UJOUJO j 

J-o 

(14b) 

T 

(14c) 

0 zz — ~r 

k 

(TuJUJq 0- 

(14d) 

(T zuj — & zujq — 0 ■ 

(14e) 


In Fig. [2l we have reproduced the values of the most 
relevant covariances as a function of k or I 3 in a case 
with T j/j Tq. The coupling k produces a shift of the 
two “temperatures” Io^uj and Ioa UoUo . It also produces 
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FIG. 3. The three eigenvalues of matrix A as a function of k 
or Jo. The other parameters are fixed as in case C of Tabled] It 
is interesting to see that the experimental observations in [ll| 
occurred in the region beyond the bifurcation, i.e. the one 
with complex eigenvalues. 


with 


B 2 


Bx 


B 0 


7o I + 7-fo 
Ho 

7o7 + k(I + Io) 


k 


Ho 

7o +7 

IIo 


(16a) 

(16b) 

(16c) 


For the solutions to be complex or not, the value of the 
discriminant must be compared to zero: 

A = 18B 2 BiB 0 - 4:BlB 0 + BlB\ - 4 B\ - 27 B$. (17) 


If A < 0 there are one real and two complex conjugated 
eigenvalues. If A > 0 the eigenvalues are all real. There¬ 
fore a bifurcation point appears when A = 0 and this 
happen depending on the values of 70 , 7 , k,I,Io. We do 
not intend to exhaust all the possibilities. In Figure |3] we 
have reported a plot of the eigenvalues when k or Iq are 
varied, keeping constant all the other parameters. The 
bifurcation when k is increased is quite evident. It is in¬ 
teresting to notice that a comparison with experimental 
data, see EH where k/I ~ 4 • 10 3 , suggests values of the 
parameter in the phase with complex eigenvalues. 

It should be noted that, even when all eigenvalues are 
real, the correlation functions (discussed in details be¬ 
low) can show non-monotonic behavior because of the 
superposition of exponentials with different characteris¬ 
tic times and with positive/negative coefficients. The ap¬ 
pearance of complex eigenvalues is a source of persistent 
oscillatory behavior in the correlation function, which is 
of course always damped at large time by the exponential 
with remaining real eigenvalue. 


the appearance of cross correlations a UUo : it is crucial to 
note, from Eg. 1141 that such cross correlation is empty in 
the equilibrium case, even in the presence of the coupling 
k > 0 . 


V. DYNAMICS 


A. Characteristic times 

The dynamics is characterized by three relaxation 
times, which are the inverse of the eigenvalues of matrix 
A. The latter are the zeroes of the characteristic (cu¬ 
bic) polynomial: a compact expression is not available in 
general. The equation for the eigenvalues reads 


- X 3 + B 2 X 2 -BxX + Bq = 0, (15) 


B. Velocity power density spectrum 


In this section we study the vpds for the probe’s an¬ 
gular velocity ui(t), which is defined as 


S(/)= lim 


1 


tror^oo 27TtTOT 


ntTOT 

/ w(f)e 27r/l df 

Jo 


(18) 


which is also equivalent to the Fourier transform of the 
autocorrelation function in the steady state (w(f)w(O)) 
(see next Section for a discussion). In the following, to 
avoid confusion with the_probe’s angular velocity w(t), 
we will use the symbol / = 2nf to denote the angu¬ 
lar frequency associated with /, whose use makes more 
compact the formula. 

By time-Fourier-transforming Eqs. (]3|) and taking the 
squared modulus, the following formula for the vpds is 
obtained: 
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S(f) 


l_ _ 7 T[k 2 + / 2 (7o ~ 2fc/p) + / 4 /g] + 7oT 0 fc 2 _ 

k 2 (7 + 7o) 2 + 7 2 {fc[(/ + /o) 2 ^ - 2 /o 7 2 ] + (7 2 - 2rt)7o 2 } + IVoP - 2 II 2 k + / 2 ( 7o 2 - 2/„fc)] + f 6 Pp ' 

(19) 



FIG. 4. Velocity power density spectrum (vpds) for some 
cases. The parameters of the model are illustrated in Table Q] 
The usual Lorentzian vpds is observed for the uncoupled case 
A, while for the coupling case with zero or small cage inertia 
(cases B and C) we recover the same shape (plateau followed 
by a bump) as in the “overdamped model” discussed in Ba¬ 
in cases D and E, that is with non negligible cage inertia 
Io, a region on the left of the resonant bump appears with 
a power law decay resembling experimental observations at 
small / [TJ . 


The formula is quite rich and may correspond to very dif¬ 
ferent shapes depending on the choices of the many pa¬ 
rameters. We are motivated by a qualitative comparison 
with the experimental shapes, see Fig. [TJ3, and for this 
reason restrict our study to a few paradigmatic choices 
of parameters listed in Tabic Q] 

In Fig. [f] the numerical computation of formula (fl9l) 
is displayed for such choices. As expected, case A (no 
coupling, typical of dilute gases) reproduces the classi¬ 
cal Lorentzian form for the spectrum of the Ornstein- 
Uhlenbeck process, as verified by substituting k = 0 in 
Eq. (fl9l) which gives 


S(f) 


Tl 

n(f 2 I 2 +P)' 


( 20 ) 


Case B is similar to the model used in (llj| : it shows 
the resonant bump due to the cage effect, but is forced 
to a plateau in the region at low /, because of the van¬ 
ishing cage’s inertia, i.e. Iq = 0 (no superdiffusion). This 
choice corresponds to an “overdamped” dynamics for the 
collective degree of freedom 9 0 , and analytically gives 


S(f) 


1 _ k 2 (Tj + Tq7o) + / 2 Ty 7 2 _ 

t r fc 2 (y + 70) 2 + f 2 [I 2 k 2 + 7 q(7 2 - 2 Ik)\ + / 4 / 2 7 q 

( 21 ) 


FIG. 5. Here, we show a comparison between the theory, 
Eq. (1191) (dashed lines) and the experimental data (symbols) 
of the velocity power density spectrum (vpds). The param¬ 
eters of the theory are T — 0.3, 7 = 1, 7 = 100, k = 6500, 
To = 70, Io = 10 4 , 70 = 10, which are close to case E. The 
experimental data come from Ref. , in the cold liquid case. 


The expresion is the analogue of the one reported in Q 
with the (not crucial) novelty that here dpt) is also sub¬ 
ject to the reciprocal of the coupling elastic force. 

Case C is similar to B, since 7o, even if finite, is still 
small. Case D and E, on the contrary, exhibit the effect of 
growing cage inertia and therefore an increasing behavior 
as / —> 0. In both cases, necessarily, S(f) becomes flat at 
very small /, since all characteristic times - even if large 
- are finite. From Fig. [4] it is clear that the effect of Iq is 
crucial: it rules the larger relaxation time and determines 
the duration of “anomalous” part of the spectrum. 

It is interesting to realize that the ramp at small val¬ 
ues of / in cases D and E has a ~ f~ 2 behavior which 
resembles the high / region. This is reasonable as co, at 
small frequencies, is somehow enslaved by the dynam¬ 
ics of u>o: the latter is, basically, a Ornstein-Uhlenbeck 
process realized at much longer timescales. This is also 
the elementary reason why the mean squared displace¬ 
ment, as discussed below and as observed in some exper¬ 
iments (nj . is close to ballistic. 

In Fig. [S] we can observe that our model - even if repro¬ 
ducing, in cases with large cage inertia (D and E), the 
ramp at small / associated with super-diffusion - does 
not do a perfect job in reproducing the whole shape of 
the experimental vpdf. Indeed, the experimental obser¬ 
vations suggest that the crossover from the cage bump 
to the ramp at small / is much smoother, basically flat: 
our model, on the contrary, predicts (independently from 
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FIG. 6. Autocorrelations C uut (t) for some choices of pa¬ 
rameters, as illustrated in Tabled] The inset shows the same 
autocorrelations at larger time, in order to highlight the per¬ 
sistent behavior of cases with high cage inertia Jo (D and F). 
Interestingly, the cage inertia lo does not affect in a relevant 
way the behavior of the autocorrelation at time smaller than 
0.1. Case E, very close to the F one (different temperatures 
do not affect the qualitative shape of the autocorrelation), is 
not shown for the sake of clarity 


the parameters) a deep elbow. It is no possible to find 
any region in the parameter space where a plateau ap¬ 
pears as in the experiment, namely, there is always a 
sharp minimum for the S(f) near the local maximum. 
Our conjecture is that, in order to reproduce the smooth 
crossover, the noise rjo(t) should be modified, possibly 
replaced by some coloured stochastic process. Notwith¬ 
standing this discrepancy, we will see in the next section 
that our theory fits very well the behavior of the mean 
square displacement. 


C. Autocorrelations and mean squared 
displacement 


all other cases with k ^ 0. The cage’s inertia Iq does not 
change in a crucial way the behavior of the first part of 
the autocorrelation, but determines a larger and larger 
persistency of correlations at late times, as seen in the 
inset. It is important to realize that the value of the 
autocorrelation at large time, even if not vanishing, is 
very small with respect to its order of magnitude at small 
times. 

From the knowledge of C uu (t), it is possible to com¬ 
pute the time-evolution of the mean squared displace¬ 
ment (msd): 


<[A0(t)] 2 ) 



dt'Xt'Mt")) 



t')C uu (t'). 


(23) 


In Figure[7]we show, for a few choices of the parameters, 
the behavior of the msd. The uncoupled case A repro¬ 
duces the standard ballistic-diffusive dynamics which is 
typical of diffusion in dilute gases. Coupling (k > 0) in¬ 
duces a dynamical arrest in the form of a plateau in the 
msd, which - later - is overcome by the slow cage dynam¬ 
ics. Such dynamics is purely diffusive in the case without 
inertia (To = 0 or small) while it is super-diffusive when 
Iq is large. When the observation time is larger than the 
time dictated by Iq (basically Iq/jo) the msd comes back 
to normal diffusion. In experiments such a very late stage 
was observed on timescale of the order of hours, see Sup¬ 
plemental Material in M- We wish to underline that, 
how anticipated, the equilibrium case F - which differs 
from case E just in the fact that the temperatures are 
equal (T = To) - displays superdiffusion as well. The 
long memory is induced by the large cage inertia Iq and 
is not related to the system being at equilibrium or out 
of it. 

To conclude this Section, in Fig.[8]we show an excellent 
agreement for the msd between theory and experiment in 
a cold liquid case. The superposition is fair along all the 
timescales. 


The two-time covariance matrix in steady state, which 
is time-translational invariance, reads 

Cij(t) d = (wi(t)wj( 0)) = exp(-At)Cij(0). (22) 

A compact expression of exp(— At) is out of question, 
as already discussed for the eigenvalues. We recall that 

d e f 

C uu (t) = C 22 (t) = (w(f)w(O)) can also be obtained by 
the inverse Fourier transform of the vpds. This is true 
for each component of Cij(t) which is related - through 
transforms - to the cross-spectrum. 

Autocorrelations for some choices of the parameters 
are numerically computed (from analytical formula) and 
shown in Fig. [G] It is evident the passage from a simple 
exponential decay of uncoupled case A (k = 0) to the 
“back-scattering” behavior, typical of cages in liquid, in 


D. Persistency. 

An intriguing counterpart of super-diffusion is a long 
memory in the velocity variable ui(t). At times shorter 
than the typical time needed to overcome the cage effect, 
the dynamics of w(t) is dominated by rapid intracage os¬ 
cillations and fast thermal relaxation. As seen in Fig. [Gl 
when cage inertia is large the autocorrelation of uj(t) dis¬ 
plays a very small long-lasting drift. Such a weak sig¬ 
nal in the autocorrelation disappears, in experiments, 
because of noise. A better observable [Ill, is the dis¬ 
tribution of inversion times for the slow part of the sig¬ 
nal, which is strictly connected to such a long memory. 
In order to make contact with the experimental results 
in (Tl| . we have studied this distribution in our theo¬ 
retical model. To remove the effect of rapid relaxation, 
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we apply a filter and study ui s (t) = i j d+T dsui(s), us¬ 
ing values of r = 2 larger than the cage relaxation time. 
An inversion time ti nv is the time between two zeroes of 
w s (f). The inversion times, therefore, represent the du¬ 
ration of “persistency” of the direction of motion. We 
remark that a theoretical computation of the distribu¬ 
tion of persistency times is a tough task which we have 
not pursued here. On the contrary, we have simulated 
numerically the model © and we have computed the 



FIG. 7. Mean square displacement for some choices of the pa¬ 
rameters, as illustrated in Table [T] The phenomena observed 
in the vpds are also observed in the MSD. After a common 
ballistic regimes, all cases with coupling (i.e. excluding case 
A) show an intermediate cage effect visible as brief plateau. 
After the plateau the cases with zero or small cage inertia Io 
(B and C) show normal diffusion, while the cases with large 
cage inertia (D and E) present ballistic super-diffusion. The 
super-diffusive behavior terminates at a time which is larger 
as Io increases. 
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FIG. 9. Distributions of inversion times for some of the cases 
in Table [T] All distributions present a - more or less visible - 
cut-off at a finite time. Looking at the coupled case (B-E), it 
is seen that when the cage inertia Io is increased, the cut-off 
time increases too. In the cases with largest cage inertia (D 
and E), the cut-off time is large enough to “uncover” some 
power-law decay. In order to present more clearly the results, 
and given that cases E and F are very close, we show only 
case E. 


distribution P(ti nv ) : which is shown in Fig. [9] From this 
plot one immediately appreciates the role of cage iner¬ 
tia Iq in inducing large inversion times. Uncoupled case 
(A) as well as the coupled cases with zero or low inertia 
(B and C) show a distribution of ti nv with a fast de¬ 
cay at a small cut-off time. When the cage inertia Iq 
is increased, the distribution develops a large tail ~ t~^ 
whose cut-off increases, while a < 2. It is interesting 
to notice that similar exponents have been observed in 
the experimental distribution, where they were in perfect 
agreement with superdiffusion predicted by a continuous 
time random walk model for ui s (t) EH- 


E. Linear response 

The aim of this section is to study the response of the 
variable ui(t) to the application of a force (torque) upon 
it. As the model is linear, the response is always linear. 
However we do not believe that the model can reproduce 
experimental observations when a large force is applied, 
i.e. in the non-linear regime it should be properly modi¬ 
fied. 

The dynamical response matrix simply reads R(t) = 
exp(-At). Therefore R uu {t) d = R^t) = It is in¬ 

teresting to notice that the same expression is obtained 
by applying the Generalized Fluctuation-Response Rela¬ 
tion discussed in 0, [TT(f: 

Sio(t) _ / dlnP(z,tu,uj o) 

MO) “ \ du 


FIG. 8. Comparison between the theory (dashed lines) and 
the experimental data (symbols) for the mean square displace¬ 
ment for the same case of Fig. [u] 


t =0 


(24) 
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FIG. 10. Deviation from equilibrium measured as a distance 
from the Einstein relation, in some cases from the list in Ta¬ 
bled] There is no deviation in the uncoupled case (A) and in 
the coupled case when T = To (F). All other cases show a 
deviation from the Einstein relation of similar order of mag¬ 
nitude, with the exception of the large cage inertia case (E) 
which is quite smaller. 


where P(z,u!,uo) is the umperturbed steady state dis¬ 
tribution, Eq. © in our case. Applying the results of 
Sec. [TV] we get, therefore: 

= Puw{u{t)u(0))+P UUo (uj(t)u]o(0))+P zu (uj(t)z{0)), 

(25) 

where we remind that /3 = a 1 is the inverse covariance 
matrix. Such an expression makes clear the role of static 
correlations, which is the crucial ingredient mod ifying t he 
equilibrium Fluctuation-Response relation [23L13714391 ] 
Indeed, from expression (12511 and Eqs. da, it appears 
that at equilibrium, i.e. when T = Tq (with or without 
coupling) one has Ruu{t) = C uul {t)/C bJbJ { 0), which is a 
way to express the so-called Einstein relation which gives 
mobility as proportional to diffusivity. Plots of the differ¬ 
ence between the response R uu (t) and the rescaled auto¬ 
correlation C uu (t)/C uu (0) are shown in Fig. [TO] When 
T ^ To (and k > 0) the response is always different 
from the rescaled autocorrelation and the Einstein rela¬ 
tion is always violated. An experiment comparing linear 
response to autocorrelation would be able to put in evi¬ 
dence the distance from equilibrium in the system. 


pared with granular gases or with slowly tapped/sheared 
(or even static) granular “solids”. The situation could 
be considered similar to the theory of liquids in the 
60’s, where neutron scattering spectra were explained by 
proposing super-simplified models [25], [HI . 

In this spirit, here we have proposed a simple model 
which reproduces two striking features of a recent exper¬ 
iment [ill ]. namely: a transient dynamical arrest (cage 
effect), and late time - almost ballistic - superdiffusion. 
A qualitative comparison with the experimental results is 
fair for the velocity power density spectrum, as it includes 
the “resonant” bump associated with the cage effect and 
the strong enhancement at low frequencies (/ —» 0) which 
seems related to superdiffusion. The comparison for the 
mean squared displacement is even more striking, and 
shows how the superdiffusive behavior emerges only when 
the “cage inertia” Iq is large. In close analogy with the 
experiments, also the distribution of persistency times is 
strongly affected by the value of Iq which controls the 
cut-off time of the distribution and the possibility to ob¬ 
serve slow (power-law) decays. We have also discussed 
how the lack of thermal equilibration between u>(t) and 
the collective variable tuo(i), which is a common feature 
of granular fluids due to non-conservative interactions, 
could be demonstrated in future experiments through an 
analysis of linear response. 

As a concluding remark, we wish to underline that the 
model proposed is purely phenomenological and lacks a 
general derivation from first principles. In particular, the 
experimental parameters such as the packing fraction or 
the average energy input do not enter in this simplified 
picture and there is no way to to deduce or estimate a rea¬ 
sonable value for Iq . The experimental observations sug¬ 
gest that, when the “granular temperature” is decreased, 
larger and larger values of Iq are realised. Tentative fits 
of experimental vpds at the smallest granular temper¬ 
atures, with the formula derived here, give the values 
of the “collective” momentum of inertia of Iq ~ 10 4 1: 
if we assume that this is the momentum of inertia of a 
solid made of a material density of steel (as the granular 
spheres) reduced by the packing fraction (order ~ 50%) 
it gives a radius of the order of 10 2 mm which matches the 
order of magnitude of the size of the container filled by 
grains in the experiment. A mechanism explaining the 
building up of cage inertia is still lacking and certainly 
deserves future investigation. 
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